Anisotropic diamagnetic response in type-II superconductors with gap and Fermi 

surface anisotropics 
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Effects of anisotropic gap structures on a diamagnetic response are investigated in order to demon- 
strate that the field-angle-resolved magnetization {Ml{x)) measurement can be used as a spectro- 
scopic method to detect gap structures. Our microscopic calculation based on the quasiclassical 
Eilenberger formalism reveals that Ml{x) in ^ superconductor with four-fold gap displays a four- 
fold oscillation reflecting the gap and Fermi surface (FS) anisotropics, and the sign of this oscillation 
changes at a field between Hci and Hc2- As a prototype of unconventional superconductors, mag- 
netization data for borocarbides are also discussed. 
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Precise determination of the node position or gap 
structure is of fundamental importance in superconduc- 
tivity study in general, especially for ever-growing so- 
called unconventional superconductors since it is indis- 
pensable in identifying the pairing mechanism for a mate- 
rial of interest. There are only a few established methods 
for the precise determination of gap structures; Angle- 
resolved specific heat and thermal conductivity mea- 
surements are notable ones Q, S S S 0- Neverthe- 
less, to reinforce the conclusion, more such spectroscopic 
experiments based on bulk quantities are desirable jg. 
As Takanaka suggested within the Ginzburg-Landau 
(GL) theory, a diamagnetic response is a strong candi- 
date if combined with an analysis of the basal plane mag- 
netization anisotropics. Needless to say a diamagnetic 
response from a superconductor is a hallmark of rigidity 
of the macroscopic wave function, containing a wealth 
of microscopic information, and it is a routine work to 
measure magnetizations to check if a material of inter- 
est is a superconductor or not. Since we are realizing 5] 
that field-angle-dependences of various physical quanti- 
ties such as specific heat or thermal conductivity should 
reflect low-lying quasi-particles around the vortex core, 
it is also expected that the magnetization contains the 
same kind of information. 

Among a vast amount of type-II superconductors non- 
magnetic borocarbides i?Ni2B2C (i?=Lu,Y) are consid- 
ered to be typical examples of unconventional ones in 
the following sense: (i) Lots of experiments have demon- 
strated the existence of gap nodes (or gap minima) 
P,I3,IM0>0|> (ii) the normal and mixed state are not ex- 
posed to the strong fluctuations such as in cuprates. (iii) 
there exist detailed magnetization measurements |llLll2j . 
To establish a spectroscopic method based on the magne- 
tization measurements it is thus important to understand 
the anisotropic diamagnetic response in these materials 
and clarify its relation to the gap structure. 

In the works by Civale et al. 0| and Kogan et al. , 
basal plane magnetizations were measured as a function 
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FIG. 1: The coordinates {x, y, z) and the crystal axes (a, fe, c) 
The induction B is rotated from the a-axis by an angle x- 



of the angle x (see Fig. between the applied field 
and the crystal axis, and it was found that the four-fold 
oscillation of the magnetization showed a sign reversal 
with decreasing the field (or temperature). Kogan et 
al. demonstrated that these behaviors can be repro- 
duced within a nonlocal London theory without quoting 
anisotropy effects of the gap function. From a theoreti- 
cal viewpoint it is quite reasonable to consider that the 
lower field (or lower temperature) behavior is well de- 
scribed within this framework since the London theory 
is believed to be appropriate for studying a vortex state 
in a field region H <C i?c2- With increasing the field, 
however, the vortices start to overlap and the physics of 
vortex cores begins to play an essential role, whereas the 
core physics cannot be captured by any London theory 
without introducing some cutoff function by hand. As a 
consequence, at least when discussing anisotropy effects 
in a vortex state, the validity of London description in 
high fields is quite unclear. If we aim to clarify whether 
the observed phenomena are generic ones or not, we need 
a theoretical approach which can correctly describe the 
anisotropy of Hc2 and the core effects. 

The purpose of this Letter is two-fold. The first 
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one is to clarify the effects of gap structures on the 
anisotropic diamagnetic response based on the quasiclas- 
sical Eilenberger formalism. The second one is to apply 
our analysis to the prototype materials, i.e., nonmag- 
netic borocarbides, and present a microscopic descrip- 
tion of the observed mysterious sign reversal of Ml {xY 
oscillation [lllll^ . Our microscopic treatment covers, of 
course, both GL theory and London theory, since the for- 
mer is derived from the Eilenberger formalism through an 
expansion about the pair field, while the latter by using 
a phase only (London) approximation. Indeed our nu- 
merical solution at finite temperatures well reproduces 
both GL behavior (M cx Hc2 — B) near Hc2 and London 
behavior (M oc \n.{Hc2l B)) in lower fields. 

Let us first explain our numerical procedure. We start 
with the following Eilenberger equation (k-Q = h = 1) 

/(e„, p, r) = (2£„ + iv • n)"^ (2.g(e„, p, r) WpA(r)^ , 

(1) 

where /(£n,p, r), /t(£„,p,r) = /*(£ „, -p,r), and 
5(£„,p,r) = Vl - /(£„,p,r)/t(£„,p,r) are the Eilen- 
berger's Green's functions. Here £„ — 2TrT{n + 1/2) 
is a fermionic Matsubara frequency, v is a Fermi veloc- 
ity, n — — iV -I- (27r/$o)A is a gauge invariant gradient, 
Tc is a transition temperature at a zero-field, and $o 
is the fiux quantum. The gap function is expressed as 
Ap(r) — WpA(r) where Wp is the pairing function with 
relative momentum p of the Cooper pair, and A(r) is 
the order parameter with center of mass coordinate r. 
Throughout this paper we treat extreme type-II super- 
conductors with large GL parameter k 3> 1 (in borocar- 
bide superconductors k ;J 10), so that the vector poten- 
tial is approximated by A = —Bzy where B = _Bx is the 
induction (see Fig. . This is indeed a good approxima- 
tion in high-K case since the correction term is of order 
0(1/k^). To solve Eq. (0), we adopt an approximation 
similar to that used by Pesch |l3 |. 
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The physics behind this approximation is that the spa- 
tial variation of / related to the phase modulation of A 
is much larger than the spatial variation of g describing 
the amplitude fluctuation. It is worth noting that we do 
not replace g in the above equation by its spatial average 
as Pesch done, in order to ensure that the correct expres- 
sion for the nonlocal GL free energy F/V given in Ref. 
Hi! are reproduced up to the quartic term. Besides the 
above mentioned justification near Hc2^ the applicability 
to lower fields is improved by requiring the sclfconsis- 
tency among /, p and g. Furthermore our scheme can 
be valid in an anisotropic case by including the contribu- 
tion of higher Landau level components to A: 
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where Aq - I.IGAT,, = V<^o/2^S, = 

{vtb/ \/TT)e~^'^'^™' , Hn is the A^-th Hermite polynomial, 
and the lengths are measured in units of r^. The real 
constants C and specify the configuration of a vortex 
lattice. Since the difference of a vortex lattice configu- 
ration is considered to be irrelevant to the quantity in 
question, we set in this paper <^ = 1/2 and v — \J a/Stt; 
the value for a triangular lattice. Substituting the above 
expression into Eq. and using a parameter represen- 
tation 
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we have 
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Here the expression for aN = e ^''^'^iI>n is given by 
Hn(z + vm — Rc A) 
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where A — {vz+ivy)p/rB- At long last we have a solution 
for /, on condition that we have the correct {d^vl-values. 
To determine {cZat} we use the following selfconsistent 
equation projected onto each Landau level: 

( H^) + 2nTj2s-')dN = ^ E ^^K/)' (8) 

^ n>0 " n>0 

where the overbar denotes the spatial average, and the FS 
average (• • •) satisfies a normalization condition (1) — 1; 
the definition is given by Eq. I|ll|) below. Our numerical 
procedure is as follows: Input initial values for {d^}, 
/, and g. Next use Eq. ^ to obtain the new / 
(and /t, g). Then iterate Eq. ((HJ to renew the {dN}- 
values, and return to Eq. ©. In order to check the 
reliability of our numerical procedure, we initially treated 
a two-dimensional case and calculated field dependences 
of each {d^} at T/Tc — 0.5. The obtained result was 
quite similar to the previous work based on the Landau 
level expansion of the GL equation (Figure 2(b) of Ref. 
0). In the following calculation we use A^ax — 6. 

The magnetization 47rM — B — H is obtained from the 
relation H = 47rVB(^"/V^). As for the longitudinal com- 
ponent Midi H), which we focus on in this paper, Klein 
et al. [T^ obtained a more convenient formula extending 
the virial theorem derived by Doria et al. 18]: 
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FIG. 2: Logarithmic field dependence of longitudinal mag- 
netization Ml at T/Tc = 0.6 for B || node (filled circles) and 
B 11 antinode (open circles). The used anisotropy parameters 
are Q = 1 and /3 = 0. Horb = 1.037$o/27rfo ilo = vo/2tvTc) 
is the orbital limiting field in the isotropic case. Inset: Field 
dependence of M^" = {Ml{x = 0) + Ml{x = 7^/4) )/2 for the 
same parameters. 



where in the above equation we approximately set 
B • H ~ Bi/ as in Ref. 0. 

Now let us assume an isotropic FS to focus only on 
the role of the gap anisotropy. In this case the Fermi 
momentum can be written as p = mvof, where f — 
(sin 6* cos sin sin 0, cos 0), vq is the Fermi velocity in 
the isotropic case, and m is the effective mass of the 
quasi-particle. As a model for an anisotropic gap func- 
tion, Wp = yj\ — acos4(0-|- x)/(l ~ acos4((?!) -|- x)) is 
used where x is the field-angle measured from the crystal 
a-axis, and is the azimuthal angle measured from the 
X-axis. Thus a denotes the degree of the gap anisotropy 
with a — \ being the nodal case. Figure El shows the 
field dependence of the longitudinal magnetization 
for B II node (filled circles) and B || antinode (open cir- 
cles) at r = O.GTc- From the inset both the GL and the 
London behaviors are clearly seen. Although the differ- 
ence of Ml between the two field-orientations is rather 
small, we can find that the at B ~ Hci is larger for 
B II node, but with lowering the field this tendency is 
reversed at a field B* . In Fig. Of a), the corresponding 
field-angle dependences of Ml(x) are plotted for several 
inductions. It should be emphasized here that the gap 
anisotropy alone can cause a sign reversal of the Mi(x)- 
oscillation. The result with isotropic FS can be summa- 
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Next we discuss the magnetization experiments [lllll^ 
for borocarbides. The observed data are incompatible 
with the above conclusion once we recall the experimen- 
tal suggestion that the nodes exist along [100] and [010] 
directions 1, The discrepancy is considered to stem 
from the unusually large FS anisotropy possessing partial 
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FIG. 3: Field-angle dependences of Ml{x) for (a) a = 1, (5 = 
and (b) a = 1, /3 = 0.4 at T/Tc — 0.6 for several inductions. 
The data at different B are vertically shifted. 

nesting [3| in these materials. As a model having such 
an anisotropic FS, we introduce a four- fold anisotropic 
dispersion 



(10) 

where a{9, (j)) = ^1-1-/3 sin^ 9 cos 4((/) + x) • Thus (3 de- 
notes the degree of the FS anisotropy. The resultant 
Fermi momentum can be written as p = {mvQ/a)f. 
The Fermi velocity v — VpCp can be expressed as v = 
Vrr + V09 + v^(f), where 6 = (cos cos 0, cos sin 0, — cos 9) 
and (j) — (~ sin(/),cos(/), 0). Here each component of v is 
given by Vr — vqct, vg = 2vq{[3/(t) sin^ 6* cos 6* cos 4(0-1- x), 
and U0 — — 2uo(/?/o-) sin^ 6'sin4(0 -I- x)- Finally, the 
area element dS of the FS divided by |v| is given by 
(iS'/|v| = {rn?VQ/(j'^)d{cofi9)d4i. Then our definition of 
the FS average is given by 
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A band structure calculation for LuNi2B2C suggests a 
rough estimate /3 ~ 0.4 so as to reproduce the ratio 
{vt)l{vlvl) = 0.128 20;] within our model. 

Starting from the isotropic FS case (/3 = 0) and in- 
creasing the /3- value, the oscillation pattern seen in Fig. 
ini^a) first tends to diminish, and when the /3-value ex- 
ceeds about 0.1, the sign of the Mi(x)-oscillation pat- 
tern is completely reversed. This is shown in Fig. ISfb), 
and the oscillation behavior well coincides with the re- 
sults of Refs. and Worth noting is that a/3 > 
in our model corresponds to the competing anisotropy 
case in the sense of Ref. 0, where the observed con- 
figuration 01 of the vortex lattice in B || c is prop- 
erly explained by the competition between gap and FS 
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FIG. 4: Logarithmic field dependence of 5Ml = Ml(x ~ 
7r/4) - MLix = 0) for a = = 0.4 at T/Tc = 0.6. Inset: 
The sign reversal field _B* in the B-T phase diagram. The 
Hc2 is determined through linear extrapolations of Ml- 

anisotropies. Note also that the main conclusion here is 
not changed by the nodal topology and effective dimen- 
sionality of a material, though the oscillation amplitude is 
quantitatively enhanced for a quasi two-dimensional ma- 
terial. 

We show in Fig. 0] the field dependence of the mag- 
netization oscillation amplitude SMl = Ml{x = 7r/4) — 
Ml{x = 0) at T/Tc = 0.6, which is to be compared with 
the experiments (Fig. 3 of Ref. and Fig. 5 of Ref. 

The characteristic behavior that 6 Ml is a slowly in- 
creasing function of the field except for a peak structure 
slightly below Hc2 is well reproduced. Finally the inset of 
Fig. 0] shows the oscillation sign reversal field B* in the 
B-T phase diagram for a = 1 and /? = 0.4. The B*{T) 
is a decreasing function of the temperature, and this is 
consistent with the experimental finding of Ref. [U The 
B*{T) could be a diagnostic quantity to characterize the 
gap function and the FS anisotropy. 

The physics behind the phenomena is explained as 
follows: Near Hc2 where the GL theory can be ap- 
plied, the anisotropy of the longitudinal magnetization 
AttMl ~ (B - Hc2)/2.32k'^ comes from that of Hc2- In 
lower fields where the London theory is more appropri- 
ate, the anisotropy of 47rML ~ (-^q/SttA^) In(7?c2/-B) 
is effectively attributed to that of the penetration depth 
A. Extending a knowledge on the relation between the 

and A-anisotropies based on the anisotropic GL equa- 
tion [i^, a naive relation i?c2(x = 0)/iJc2(x = 7r/4) ~ 
'^^(x = 0)/'^^(x = ^/4) is expected to hold, where A(x) 
means an effective penetration depth perpendicular to 
Hc2{x)- Namely if Hc2 is larger, then 1/A^ is expected 
to be smaller, and this causes the sign reversal of the 
Ml (x)-oscillation. 

Before ending we briefly discuss the origin of the pecu- 
liar angular variation of Ml(\) seen in Refs. andll^ 
The low field Ml (x) [ill UM have sharp maxima around 



[100] and [010] directions, and broader minima around 
[110] and [110] directions. The sharp maxima remind us 
of the observed cusp-like minima of thermal conductiv- 
ity and specific heat jllQ around [100] and [010] direc- 
tions. These have been argued as a result of so-called 
"s+g" pairing function Q. We calculated Ml{x) for the 
pairing Wp oc (1 — sin* 0cos4((/() -I- x)) with point nodes 
keeping P — 0.4. However no such characteristic struc- 
ture was observed in the angular dependence of Ml{x)- 
This means that explanations for the observed structure 
of Ml{x) may need a different mechanism. 

In conclusion we have microscopically studied field- 
angle-resolved basal plane magnetization oscillations, 
and demonstrated not only that a careful measurement 
of the magnetization can be a potentially useful tool to 
identify the nodal position of the gap function when a 
material of interest possesses less anisotropic FS, but also 
that the experimental data for borocarbides are well re- 
produced by considering both gap and FS anisotropies. If 
combined with other field-angle-resolved quantities [j^ . 
such as specific heat and thermal conductivity at low 
temperatures, we can further narrow down the possible 
pairing symmetry or gap anisotropy in various conven- 
tional and unconventional superconductors. The lesson 
from borocarbides tells us that for materials with unusu- 
ally strong FS anisotropy such as borocarbides we should 
be careful to judge the nodal position through magnetiza- 
tion measurements, since unusually strong FS anisotropy 
can reverse the conclusion. 

After the submission of the first manuscript we learned 
about a preprint by Kusunose 113] which study the effect 
of the gap anisotropy alone based on a simplified version 
of our treatment. The result is consistent with ours when 
only gap anisotropy is considered. 

We would like to acknowledge useful discussions with 
T. Sakakibara and Y. Matsuda. 
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